Quasirelativistic electronic T-matrix for the short-range perturbation 
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Exact T-matrix for the delta-function short-range perturbation of the (3+l)-Dirac equation has 
been derived. Separability of the potential in the angular momentum representation is used. A 
characteristic equation for the T-matric poles determining the bound and resonance states has been 
obtained. The calculated T-matrix determines transport properties of the narrow-gap and zero-gap 
semiconductors . 

I. INTRODUCTION 

Transport theory and theory of the electronic spectrum need in calculation of the T-matrix including its " nonphys- 
ical" part outside the mass shell. This makes it necessary to consider exactly solvable models of perturbation. One 
of them is the delta function potential. In the case of the two-band nonrelativistic problem described by the Dirac 
equation the problem of bound and resonance states was considered with a use of such potential in the paper ??. The 
short radius potential S (r — ro) was considered. This potential has no a singularity at r = and is separable in the 
angular momentum representation. Here were consider the problem of the bound and resonance states with slightly 
different approach and calculate an exact T-matrix. We take into account possible difference of the perturbation 
matrix elements calculated on wave functions of the upper and lower bands that is equivalent to consideration of both 
potential and mass perturbations. 

II. CHARACTERISTIC EQUATION 

As it is well known, the electronic states near the band edge are determined by the quasirelativistic Dirac equation 
jllj . It can be written in the momentum representation in the form: 

s<?PX (P) + (™s 2 -E)<p(p)+f d 3 p'Vi (p - p') if (p') = 0, (1) 

s<j W (p) - (ms 2 + E) X (p)+J d 3 p'V 2 (p - p') X (p') = 0, (2) 

where a is the Pauli matrix, s is the quasirelativistic limiting velocity of band electrons, 2ms 2 — E g . The potential 
Fourier transform Vi (p) can be expanded into a series of the Legendre polynomials Pi (cos 9) in the case of the 
spherical symmetry [2j: 

Vi (|P - P'|) = V l Pl ( cos0 ) ( 2l + !) / 47r ' ( 3 ) 
i 

Vi(p,p) = 2 drrV t {r) = , (4) 

Jo VPP 

where J; + i/2 (jpr) is the Bessel function. The wave functions can be expanded into a series of spherical spinors flji m (n) 
1: 



X% m (p) J \-gE(p)i- l 'n j ,, m (n))' 



(5) 



where I = j =p 1/2, I' = j ± 1/2, j is the total angular momentum quantum number, n = p/p, 8 is the angle between 
the vectors p and p', while 

anQji m (n) = i l ~ l £ljv m (n) . (6) 



Signs ± correspond to x = ± (j + 1/2) . Thus, eqs. ((T|), ^ take the form 



2 



(ms 2 - E) f E (p) n jlm (n) - spg E (p) (-1)' l ' Sl jlm (n) + / d 3 p' £ v£ (p,p) P h (cosfl) f E (p') n iJlTOl (n') = 0, (7) 



l\ , ??1 1 



(ms 2 



E) g E (p) Sl jUm (n)+spf E (p) (-!)''-' tt^m (n)- ^ d 3 p' ^ K,' 1 (p,p') P Zl (costf) ,g B (p') fi iIimi (n') = 0. (8) 



Making use of the addition theorem for the spherical functions 



47T 

m 

of the spherical spinors orthogonality [ij 



£ ^ (n) Y lm (n') = — ^ (costf) , (9) 



y rfo^*/m ( n ) ^'liiTOi (») = 5j,jA,h S m,mi ( 10 ) 



, / j+l/2±(m-l/2) y- / s 

/ \ I ± V 2(j + l/2 T l/2) ^-1/2 W 

n J=I±1/a , lim (n) - i / j+1/2T(m+1/2) I (11) 

V 2(j + l/2 T l/2) r ''™-l/2 W 

we obtain from (??), (??) 

/>OC 

(ms 2 - E) f E (p) + spg E (p) + / dp' (p'f f E (p 1 ) V? (p, p') = 0, 

Jo 

/>00 

(ms 2 + E) f E (p) + spg E (p) - / dp' (p'f g E (p') V/ (p, p') = 0, (12) 



where i — Z' = ±1. 

Zero-radius || and separable potentials [|| are popular in the nonrelativistic scattering theory. However, the Dirac 
equation is extremely sensitive to a singularity of the potential [? ] . The singularity of the delta function potential 
can be regulated putting the delta function to the spherical shell 

Vi (r) = ^S(r - r ). (13) 

Substituting (??) into {4| we obtain the separable in the angular momentum-momentum modulus representation 
potential: 

V? (p,p')=v\(p)v\(p'), (14) 

where 



^)=y^w*-o). as) 

Equations ([12]) become degenerate and can be written as follows: 

POO 

(ms 2 - E) f E (p) + spg E (p) + v[ (p) / dp' (p'f f E (p') v[ (p') = 0, 

Jo 

roc 

(ms 2 + E) f E (p) + spg E (p) - v\ (p) / dp' (p'f g E (p 1 ) v l 2 (p') = 0. (16) 

Jo 

Introducing the functions 

/>oo />oo 

F(E)= dpp 2 f E (p)v[(p), G(E) = dpp 2 g E (p)v l 2 (p), R(p) = (s 2 p 2 + m 2 S 4 - E 2 y\ 
Jo Jo 
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we obtain the algebraic equation set 



F = G dpp 2 sR (p) v{ (p) v l 2 (p) -(E + ms 2 ) F / dpp 2 R (p) (v[ (p)) , 
Jo Jo 

G = F ^ dpp 2 Sj R (p) «i (p) 4 (p) + - ms 2 ) G J dpp 2 R (p) [v l 2 (p)J 
The solvability condition for this equation set gives the characteristic equation: 

1 



, E + ms 2 Vi , , , , 

1 H : 4fl/2 (wo) Kl+i/2 (Kr ) 



h 2 s 2 Anro 



E-ms 2 V 2 ° _ . . . , 



ft 2 s 2 47rro 



167T 2 Vq?1 2 S 2 



where k 2 = (m 2 s 4 — -E 2 ) fi. 2 s 2 , J n (x) , 2f n (a;) are the modified Bessel functions. 



III. CALCULATION OF THE T-MATRIX 



(17) 



(18) 



Let us determine the standard representation bispinor basis as follows 



(r,p|p,X, + ) 



y/ Ep + mR pi (r) (r/r) 
V -y/Ep-mRpv (r) n v (r/r) 



(19) 



/ I _ s _ » / y /gp - m -Rpj (r) (r/r) 
1 j \,^T^Hp,,(r)n^(r/r) 



(20) 



where A = ±, x = (j, Z, m); Ep = yj m 2 s 4 + p 2 s 2 ; I — V = sgn x; x = ±1, ±2, . . . ; R p i (r) = (P r ) > 

g is the bispinor index. We work in units with fi = s = 1. The exact wave function of the out-basis satisfies the 
Lippmann-Schwinger equation: 



(r,g\p, X ,\) = (r,9\p,xA)-Yl J d * rG l a ' ('"OW MWb.X.A) 



The perturbation matrix can be written in the form 

V (»•) = ^ v 



(21) 



(22) 



Distinct values for the upper and lower bands stem from different symmetries of the Kohn-Luttenger basic functions. 
We define Vi — V2 — V+; V3 = Va = V-. The T- matrix can be defined os follows: 



(p', X ',\'\T\p,x,\)=J2 J d 3 r(p\x',\g\r,g)V g (r) (r,g |p, X ,A) 



(23) 



Using (|21|) and expanding the Green function into a series of the Dirac Hamiltonian cigcnfunctions we obtain the 
equation for the T- matrix: 



(p', X ', A'| T (E) |p, X , A) = (p', x, A'| V (r) |p, x , A) - 

(p', X ',X'\V(r) \puXlAx) 



J E 



E Pl sgn\i — E — iO 



(p liX i,Xi\T(E) |p,x,A). 



The perturbation matrix elements are determined as follows 

(p',x',A'|V(r) |p, X ,A)=^ /" d 3 r(p', x ',A 5 |r, ff )V/ g (r, 5 |p, x ,A) 



(24) 



(25) 
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This matrix takes a form of a sum of the two factorized expressions in our case: 



(p', X ', A'| V (r) \p, x, A) = / do [(p', x', A'| r , n, +) (r , n, +| p, x, A)] + 
(p',x',A'| r ,n,-) (r 0) n,-|p,x,A) VI. 
We seek a solution of the Lippmann-Schwingcr equation in the form: 



f y ^+^J ^i (r) (r) 0„ (r/r) \ 
/2E; I -^E p -mR pV (r) G^ (r) (r/r) i ' 



(26) 



(27) 



(r,3 b,X,-> = 



y/W^Rpl (r)F { pl } (r) ^ (r/r) 



pi 

.(-) 



/2^p y ^ + mRpv (r) G^ (r) (r/r) 
The T-matrix can be expressed in terms of the and G ± functions as follows: 



(28) 



{p uX u+\T{E) \p, X ,+) 



5 xx' 



(E V1 + m) (E p + m)V°F^R; ih (r) R pl , (r) 



2\/F p E p 

y/(E Pl - m) (E p - m)j V°G p +?R* pill (r) R pl , (r) , 



(29) 



( Pl ,Xi,-\T(E) \p, X ,-) = 



S xx' 



2yJ E p Ep 



yJ(E Pl - m) (E p - m)V°F£R; ih (r) R pl% (r) + 



^{E P1 + m) (E p + m )V G pl ?R* pil , (r) R pl{ (r) 



(30) 



(p 1)Xu -\T(E) \p, X ,+) = 



<V (-1) 



3+1' 



(E P1 ~ m) (E p + m)VfF^R; ih (r) R ph (r) 



2\/E p Epi 

\J {Epi + m) (E p - mjy'G^j, (r) i? p/l (r) 



(31) 



{p 1 , Xl ,+\T(E)\p, X ,-)- XX ' 



(E P1 + m) (E p - mWlF^ R; ih (r) R ph (r) 



l^jEpEpi 

yj (E P1 - m) (E p + m)V°G { -)R; iVi (r) R ph (r) 



(32) 



Substituting (j2"5)) , ([3H)l , (f3"Tj) and (|3"2")) into (f2~4")l we obtain two independent sets of equations for the pairs of functions 
F(+\ G« andF(-), GH : 



~ / F {±) 



C P!P V + U P1P ' C p!p V - U PlP 

c p lP v +°PxP C pip V_0 PlP 



where 



„±± 



\J{E Pl ±m)(E p ±m), c p ^ = x/ {E Pl ± m) (E p =p m) , 



(33) 



(34) 



b p ^ ] =K lh Rph, 6fcp _) =K^ R pI'v b p\v ] =P R * Pl hM> =pR; il ,R ph . 



(35) 
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Matrix elements of A read 



V?V° 



An = V+&^ [1 + V+ ( E + "0 B++ ( E )) ± ^&p~ip )b±T ( E ) 



A ± 



V°c*%--) [1 + V* (E m) B- (25)] ± Y±^ c TT b (+-) b ±t {e) , 



A ± 

A 12 



P 

yoyo 



y° c ^6(;~) [i + y_° (is - ™) b~ (25)] ± ^- c t±6(+-) b ±t {e) 

= VZcgbgp [1 + Vl (E + m) B++ (E)] T Y±Y^^± b i-+) B T± {E) . 



where the matrix B (E) is determined as follows 

B{E) = 



dp 



E 2 -E 2 - iQ' 



Equating the matrix A determinant d to zero, we obtain the characteristic equation, which was derived in [? 
in the Section II using a different approach: 



(36) 

(37) 
and 



d (E) = det A : 



l + V?(E + m) l%L 



E 2 -E 2 - iO 



l + V°(E-m) dp 



E 2 -E 2 - iO 



E 2 -E 2 -iOJ ^ E 2 -E 2 -iO 
Solving the inhomogeneous equation set (|33[) we obtain an expression for the T-matrix: 

<5> 



(pi,Xi s Ai|T(£) |p, Xs A) = 



where the matrix elements of t are determined as follows 



'XX l 



2d^JE p E pl 



t\ ± \ (E) 



<.(+-) 



(±| 1 1±) - c ±± y>(++) [1 + F° (25 - m) B~] T c T±l/ + K_ ^i^i3 ± T + 

p 



(38) 
(39) 

(40) 

(41) 
(42) 



( ±l t |t> = (-1/(0 J cT±y +° 6 ^ +) t 1 + y -° ( s - ™) ± c++^v_°^2?*±- 

y°c ±= F [1 + Vf (E + m) B++] T c TTy0y0& B ±T. 



(43) 



The formulae (|4"0")l . (|4T|) , and (14"3")) give an exact solution for the T-matrix for arbitrary p, pi, and E. Taking these 
values on the mass shell p = p\ , E = y/m 2 + p 2 , we obtain the scattering matrix. 



Conclusion 



Non-relativistic problem of the electronic spectrum and scattering described by the Dirac equation is considered in 
the case of the two-component short-range perturbation. An exact T-matrix both on- and off-shell has been calculated. 
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